1/f Spin Noise and a Single Spin Detection with STM 
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We propose a novel mechanism for single spin detection based on the 1/f spin current noise. We 
postulate the 1/f spin noise for the tunneling current, similar to the ubiquitous 1/f noise in magnetic 
systems. Magnetic coupling between tunneling electrons and localized spin S then leads to the peak 
at Larmor frequency in the power spectrum of the electric current fluctuations / 2 . The elevated 
noise in the current spectrum will be spatially localized near the magnetic site. The difference in the 
power spectra taken at the Larmor frequency and elsewhere would reveal the peak in the spectrum. 
We argue that signal to noise ratio for this mechanism is on the order one. 

In addition we discuss the asymmetric lineshapes observed regularly with this measurement. 
We show that such lineshapes are in accordance to the random sampling done with the tunneling 
electrons. Yet, this predicts a linewidth at least one order of magnitude larger than observed 
experimentally which is likely to be due to electrostatic repulsion between the tunneling electrons 
and temporal correlations in the tunneling process. 

PACS numbers: 76.30.-v, 07.79.Cz, 75.75.+a 

1/F SPIN NOISE 

The phenomenon of 1/f (flicker) noise is known for almost 80 years 0. It describes the deviation from the flat 
spectral density expected from a current made of uncorrelated charge carriers - at low frequencies. In this range, 
the spectral density was found to obey a power law of the form 1// Q where / is the frequency, and a = 0.5 — 1.5. 
Flicker noise appears in endless number of electronic devices, in music Q in ocean streams (HQ an d in many entirely 
different systems. This is one of the most universal phenomena, yet, one of the largest enigmas in physical sciences. 

An early explanation to this phenomenon was that the 1/f noise can arise from a superposition of relaxation 
processes |fj. In this model the noise is described as a superposition of consecutive random events, each starts at a 
certain time to and follow a simple exponential relaxation law: N(t — to) = Noe~^~ to ^ T . The power spectrum of one 
such event is a Lorentzian. The power spectrum of a large number of such consecutive random events, all with the 
same r is also a Lorentzian. If, on the other hand, there is a distribution of relaxation times P(r) ~ 1/r, from n 
to T2 than the overall spectral density will obey a power law ~ 1// in the range between Tj -1 and t^ 1 . A common 
denominator to many mechanisms proposed for different 1/f phenomena 0, is a distribution of relaxation times. 

1/f noise is ubiquitous in STM tunneling current although its origins are a source of continuous mystery that is not 
fully understood until now 0, 0, 0, 0] . Unlike the well known shot noise 1/f noise is proportional to the square of 
the current such that: (I 2 (uj))/I 2 = const. Where {I 2 (to)) is the spectral density of current fluctuations (in units of 
A 2 /Hz) and I is the current. Empirically, the current fluctuations in 1/f noise are known to obey the Hooge formula 
|l2j {I 2 {uj))/I 2 = al/[fN] where N are the number of current carriers in the sample and a is of the order of 0.01. 
The appearance of 1/f noise in the STM is a surprising observation because the normal "explanation" of this noise is 
a fluctuating defect with a wide distribution of relaxation times. Such explanations are not suitable here, because of 
the extremely local nature of the measurement. The STM measurements of 1// noise (in a voltage of — 0.5U) gives a 
peak at zero frequency where the intensity of the current noise (I 2 (to)) is 20 times larger than the Thermal (Nyquist) 
noise which will give the order of magnitude of the noise in higher frequencies (at room temperature) and amounts 
to 41/ ' Aj\fWz . The width of the observed 1/f noise peak at zero frequency is of the order of 10-100KHz. 

A special case of 1/f noise relevant to the present paper is the magnetic flicker noise. In this case the noise is 
due to fluctuations in the magnetization. Normally, it is assumed that these fluctuations constitute an additional 
noise source in any magnetic system. However, low frequency magnetic noise is hard to observe experimentally and 
sensitive detection techniques and special systems are required. Such fluctuations were detected for the first time in 
spin glas ses with a SQUID magnetometer [13|,1lJ. They were detected with SQUID also in antiferromagnetic thin 
films |l5j and on superparamagnetic nanoparticles . Other ways to detect 1/f magnetic noise is by following the 
resistance fluctuations close to a certain ferromagnetic transition (colosal magnetoresistance) , and to observe the 
electrical noise generated in small Hall probe contacts 0. The 1/f noise generated in the Hall contacts is much 
larger than the usual 1/f noise measured in the same system which is unrelated to magnetic fluctuations. Despite the 
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failure to measure spatial correlations it is obvious 01 that there is a certain coherence length such that two noise 
measurements done at a distance smaller than this length will give the same results. 

Although we have no direct experimental evidence we claim that since conduction electrons in metal constitute 
also a (special) magnetic system, there should be 1/f magnetic fluctuations also in the spins in such a system which 
is also paramagnetic. ESR of conduction electrons is known for many years 0,H(|. These spectra are known for 
their special lineshape [2l| and for the narrow lines (Ti and T2 are equal and large). We would like to emphasize 
that the sequence of individual dephasing events which are responsible for T2 relaxation of conduction electrons are 
very similar to the sequence of relaxation events which are responsible for 1/f noise in general. Therefore, there are 
distributions of relaxation times that are expected to give 1/f magnetic fluctuations in conduction electrons in metals. 

Such 1/f magnetic noise is expected to be of general significance in spintronics applications. We discuss here the 
implications of such 1 /f magnetic noise (cither in the spins of the conduction electrons or in a regular paramagnetic 
system) on single preccssing spin detection with STM (ESR-STM). It will be shown below that this 1/f spin noise 
explains the observation of the signal from non magnetic tips and elucidates several aspects in this technique. We 
find that the interaction of the tunneling electrons spin with the local impurity spins simply couples the 1 /f magnetic 
noise with the noise of the local spin S thus creating a peak in the current noise at Larmor frequency. 

ESR-STM is a technique that is using the extremely local nature of the STM measurement to detect the precession 
of isolated spin centers on the surface. When a tip of an STM is located above a paramagnetic spin center (in the 
presence of an external magnetic field) the tunneling current is modulated by the precession. It was shown [22I E3 | 
that the AC current at the Larmor frequency is spatially localized within 0.5 — lnm. It is the spatial localization that 
indicates (though it must be proved) that this technique is capable of detecting a single spin. In addition it was proved 
that the frequency of the signal is dependent on real time on the size of the magnetic field fiil . More recently 
similar experiments have been done on the paramagnetic BDPA molecule l2a. A recent paper shows that ESR-STM 
can be done also on a TEMPO molecule, revealing the hyperfine spectrum [39j . The interest in this technique has 
risen sharply recently due to the possibility to manipulate and detect a single spin [23, and due to the possibility 
to use it for quantum computation |27t l28j | . There have been many proposals for the mechanism of this phenomenon 
[29L IsnL I31L l^a. Is^. Is^. I35L ^6|[ ■ However, an experimental verification for any of the proposals is still required. 

In our previous papers [3(1 l3~fl | we discussed the following question: What is the role of the Heisenberg exchange 
interaction in ESR-STM? Under which circumstances can a tip emitting tunneling electrons with a random spin 
orientation create an elevated noise level at the Larmor freuqency through interaction with the single processing spin? 

It was argued in [8fl 1.3 lj that existence of long time correlations in the temporal spin polarization of tunneling 
current is sufficient to provide the elevated noise at the Larmor frequency in the current noise. Here we further build 
upon this idea and show that the 1/f noise in the tunneling spin current is sufficient to produce the effect. Wc stress 
that no correlations in the spin polarization of the tunneling electrons within the precession period is required to 
produce the elevated noise at Larmor frequency in this approach. 

We discussed the spin dependent tunneling matrix element: 

G = Go exp[-[(F - JS(t) ■ s(t))/F ]]. (1) 

s is a spin matrix with implicit spin indices that are given by conduction electron spin operator s l = l/2a l a g,i = x,y,z 
, S(t) is the local impurity spin. In an external magnetic field it will have a random dynamics with remnant of the 
precession at the Larmor frequency. F is the barrier height (typically 4eV), J is the exchange coupling and Fq is 
the energy related to the distance d between the tip and the surface Fq = d 2 /8md 2 . In this term, the tunneling 
probability depends on the relative orientation between S and s. We write: 

G = Go exp[-[(F/F ]](cosh[(JS/2F ) + s(t) ■ n(t)sinh(JS/2F a ). (2) 

Where n is a unit vector in the direction of S. 

The Hamiltonian we consider describes a spin dependent tunnneling matrix element between the tip (L electrode) 
and the surface (R electrode) 

H = J2 e(k)cL(k)c ia (k) + (L -> R) + £ c [ Q (k)G a/3 c fl/3 (k') (3) 

k,o k,k' 

We assume that the magnetic field is along z axis: B\\z with corresponding Larmor frequency ujl = g^isB. The 
tunneling current operator will contain the spin independent part that we omit hereafter and the spin dependent 
part: 

6I(t)=G 1 n(t)I s (t), (4) 
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FIG. 1: Schematics of the ESR-STM experiment is shown. The fluctuations in the spin polarization of the tunneling electrons 
at the time scale of the precession T will be nonzero and will scale as <—. Once the tip is positioned close to the localized 
spin, the exchange interaction between the localized spin and the tunneling electrons will modulate the tunneling current. The 
conditions in which this random modulation of the spin tunneling current with 1/f spectrum will create a ujl peak are discussed 
in the text. 



where G\ = Go exp[— [(F/Fo]] sinh[^=r], and I s (i) is a spin current between tip and substrate. The dc current at a 
given bias V is : Iq = goV, go — Goexp[—[(F/Fo]]cosh[JS/2Fo\. The current-current correlator, normalized to dc 
current is then: 



= (sinh[ ^p 2 



E 



<n*(iK(t')>/j(W) (5) 

More rigorous treatment, similar to the one done in |37j |. where one takes into account the effect of the local spin 
on the mass current in the ignored terms will lead to the contribution of the same order as the term we are focusing 
on. Hence we use a simplified formula for the current that gives the right order of magnitude estimate. 

The change in the tunneling conductance due to exchange interaction between the tunneling electrons and the 
localized spin leads to Eq(4): SI(t) ~ I s (t) • S(t). Only the transverse component contribute to this term (S x (t) and 
S y (t)). When the spins of the tunneling electrons are completely uncorrelated we can say that the dispersion of the 

tunneling current which is dependent on the spin (over one precession period T) is: X^(^ 2 ) ~ Eili [ s x(~ti)Is,x(ti)] 2 + 
(x — > y) ~< N > . Or, the relative dispersion compared with the tunneling current will be: < I 2 > /I 2 ~ 1/ < N >. 
To estimate the magnitude of the spin dependent dispersion, namely the size of the noise which is due to the interaction 
with the precessing spin, we use: 

(< I 2 >) 1/2 /Io ~ 2/Nsinh[(JS/2F ). (6) 

If we take a typical value of J = O.leV, F = AeV, Fq = 0.1ey(for d = AA) and S = 1/2, we get a value of 0.01, namely 
rf intensity of 10 picoamperc, which is within the right order of magnitude. 

The spectrum is detected at the frequency domain, and by looking at the Fourier transform of the current fluc- 
tuations. This means that the spectrum is a convolution of the power spectrum of the single spin fluctuations: 
< S 2 (uj) >= r^^jWppi and the power spectrum of the spin tunneling current: < I s> i(w)I St j(— oS) >= Sij < I 2 (lu) >■ 
Here, to be specific, we assumed the spectrum of a localized spin to be Lorenzian, however the argument is applicable 
to the case of general spectrum < Sf (w) >. 



E 

i—x,y,z ' 



iifl = (sinh[— l) 2 
I 2 Q 1 [ 2F U 



If the tunneling electron spins are completely uncorrelated, then the power spectrum is white noise, and the expected 
signals will be scattered over the whole frequency range. In order to get a peak at the Larmor frequency, some 
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correlation at zero frequency is required, namely the tunneling electron spins have to have a some temporal spin 
polarization in the long time limit [35j. 

To give an example we consider the case where several magnetic atoms are adsorbed on the tip. Then the spins 
of the tunneling electrons will be under the influence of slowly changing spins h; of magnetic atoms due to exchange 
interaction Hi nt = J2i Jihi(t)(r(t) = H(i)cr(t) where Ji are the exchange couplings to the i-th atom and we define H 
as a sum of random fluctuating magnetic moments of individual magnetic impurities adsorbed on the tip. Here 
we assume that only the "field" value H at the end of the tip matters as the tunneling electrons are affected by this 
field the most before tunneling out of the tip. This fluctuating exchange field will spin polarize the tunneling current. 
The correlator is < H(t)H(t') >= E, < hi(t)hi(t') >. The correlation function from each impurity assuming they are 
independent, obeys < hi(t)hi(t') >= e - '* - * ^ ri . Average over distribution times Tj will with probability distribution 
P(r) ~ 1/r between t\ and T2, will yield the 1/f noise for the field fluctuation power spectrum. This follows directly 
from evaluating the Lorenzian in frequency space ^t^+i > avera g e< i over distribution of relaxation time r with P(r); 
e.g see Ref[7]: 

(H 2 H) = / 2 drexp(ia;(t-t , ))P(T)(h i (t)h i (t')) ~ - (8) 

Jti u 

Hence the tunneling current will acquire spin polarization fluctuations. It is natural to assume that the spin fluctua- 
tions of the tunneling electrons will be proportional to fluctuating exchange field H. The degree of spin polarization 
of the tunneling electron is proportional to the Zeeman energy (<7^bH)/W related to the total bandwidth W ~ 1/Nq. 
We get 

I a (t) = Igix B n(t)N (9) 

here we assumed that the effect is proportional to the tunneling current /. Indeed tunneling electrons will sample the 
random field H at a rate at which they tunnel and this rate is given by the electric current /. 
We now consider the specific case when spin tunneling current I s has 1/f noise component. 

{I 2 {u)) = £-C = aP (10) 

for all components i = x,y, z. The magnitude of the 1/f noise is to be given by I 2 up to an unknown numerical factor 
a |l2j . where N is absorbed in a now. 1/f noise is peaked at zero frequency and will provide a peak at the Larmor 
frequency in the convolved spectrum, Eq.(JZJ). We also remark here that in a similar fashion one can get the 1/f noise 
in the spin current due to 1/f noise in the electronic current in a presence of a constant polarizing field. We will not 
consider this possibility here. 

Spin tunneling current will have white noise asymptotic at high frequencies. At lower frequencies, similar to the 
unpolarized electric current 0, 0, Hoi 1 1 1 1 ] , the relaxation processes that control spin relaxation in the current will 
contribute. To give a physically plausible argument about the origin of the 1/f noise consider relaxation time of the 
tunneling spin current: 

(I s 4t)I SJ (t')) = Stj expHi - f\/r) (11) 

this correlation function leads to the Lorenzian for the noise spectrum with the width 7 = 1/r: (I 2 x (uj)) = 7r ^ 2 ' 7 +7 2) ■ 
Now we assume that there is a distribution of relaxation times with probability distribution P(r)dr. The possible 
origins of the distribution of relaxation times in tunneling current can be the multiple other spins in the vicinity and 
in the tip that produce spin relaxation with different times. To obtain 1/f noise in the current correlator one has to 
assume that P(t) ~ 1/r in some window of relaxation times, as was mentioned earlier. 
The power spectrum of the electric current, using Ea. (|10|) is: 



r JS ]\2 1 

l2Fo" [(w-w L ) 2 +7 2 ] 1 / 



(sinh[^]) 2 T7 — (12) 



Where we used the low frequency asymptotics for spin current Eq. (|10|l and assumed Lorenzian for spin correlation 
function. We get finally : 

(I 2 ) 1 

i^^ a/n W^yTW 2 (13) 
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Where we introduced the detuning parameter 8lo = lu — w^. 

The signal to noise ratio is controlled by the parameter a and is: 

(») 

Without the microscopic model it is impossible to know exactly what is the magnitude of a, N, e. Self consistent 
treatment would give the scattering rate 7 ~ 10 6 Hz that will be determined by the maximum of the intrinsic decay 
due to backaction or the extrinsic decay due to environment. As a guidance we take available STM data for the 
tunneling current noise H H [H El. In these data the l/f noise is clearly seen above the high frequency noise 
floor below 100kHz. From Hodge formula we estimate a ~ 10~ 2 . We take TV ~ 10 similar to mass current case, as 
estimated in [Tlj . It is convenient to relate the power spectrum of the current fluctuation at the peak at u> = cjl to 
the shot noise power spectrum (I^ hot (uj)) — el. S/N ratio is in this case is bound from above by the number on the 
order of unity 



(i: 



where we wrote explicitly the dependence of the current fluctuation on TV 



- 0(1) (15) 



RANDOM SAMPLING 



In the previous section we discussed the conditions in which the random orientation of the spins of the tunneling 
electrons can, through interaction with a single precessing spin, give a signal at the Larmor frequency. Nevertheless, 
this is only part of the picture. Since the tunneling electrons are probing a periodic precessional motion adiabatically, 
sequential tunneling and temporal correlations are essential, as will be discussed below, for getting a signal with a 
narrow lincwidth. 

One of the most important characteristic of the signals detected with ESR-STM is their lincwidth. the lincwidth 
might reflect both the lifetime of the spin state, and the back action effect of the tunneling electrons on the processing 
spin. However, it is necessary to discuss an additional source of the linewidth which is directly related to the way in 
which we probe the precession. We claim here that the periodic precession of the spin is sampled by the tunneling 
electrons. For each electron the precessing spin looks static (Adiabatic process). The simplest case which will be 
discussed here, is that the tunneling times of the tunneling electrons are uncorrelated. (As will be shown this is not 
a realistic assumption). In this case we can say that the sampling times obey an exponential distribution (Poisson 
process) . We show here using the l/f spin noise model described above that this is in accordance with the asymmetric 
lineshapes observed in the experiments. Also such a sampling process, should lead to a rapid increase in the linewidth 
when the magnetic field is increased. The data published so far, show such a trend. 

We recall that in a tunneling current of In A there are 6.25 x 10 9 electrons/sec. The measurements were performed 
in a frequency of 2 — 8 x 10 8 . This implies a ratio (Rt) of 0.033 — 0.13 of the average time between two tunneling 
electrons and the precession time. According to the Nyquist sampling theorem it is impossible to perform proper 
sampling of a periodic function if the sampling dwell time is above half of the time that it takes to complete one period. 
This is a fundamental limitation because it makes no sense to increase the field such that the precession frequency 
rises above 3 x 10 9 Hz (for a current of one nanoampere). However it is anticipated that this should be important 
also in much smaller frequencies. Due to the uncorrelated nature of the tunneling current, it is expected that for a 
periodic function that is sampled in random, (according to the Poisson distribution) increasing the frequency means 
that more and more sampling times (between two consecutive tunneling events) will be larger than half of the period. 
This implies that at larger fields (or smaller currents) we expect that the lincwidth will increase. Moreover even in 
one peak we anticipate that the high frequency side will be broader than the low frequency side, which implies an 
asymmetric lineshape that has a larger slope on the low frequency side. This is precisely the behavior that we see in 
the measurements. 

Figs. 2,3 shows two spectra that were published already. Fig. 2 shows a frequency modulated signal that is taken 
from Ref. |24j . The splitting in this signal was used to prove that the frequency depends on real time on the value 
of the magnetic field. However it is easy to see that the original spectrum is asymmetric and distorted. A similar 
distortion is seen in the (already published) spectrum (Fig. 3). This spectrum was taken from Ref. [2(|. The shape 
of the spectrum is completely asymmetric and looks very similar to the Poisson distribution function. As discussed 
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FIG. 2: A frequency modulated signal in a field of 150G. The field modulation parameters are A H = 27 mG and the modulation 
frequency is 300KHz. Two spectra are presented. The upper figure represents the data from the STM, lower is simulation. It 
is clear that the original line-shape is not symmetric. Taken from |2^| 
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FIG. 3: An asymmetric ESR-STM spectrum of a BDPA cluster in a magnetic field of 210 G. Taken from |2£ 



in details in Ref . |2,"J| , this asymmetric lincshape is the cause of getting an absorption lincshape when phase sensitive 
detection is applied. The asymmetry was explained there as a rapid passage phenomenon. As will be shown here, 
we believe now that it can be explained by the physical nature of the measurement. It is important to emphasize 
that this is not always the case. Namely other lineshapes were also observed due to spectral diffusion. However, this 
represents the majority of cases. 

In order to explain this trend we have done very simple simulation. In this simulation we have estimated the value 
represented by Eq. 8. Wc have taken a sinusoidal function for S and a random function for I s which has a long time 
correlations (1/f peak). However, in this simulation this function is sampled in times that are determined by the 
exponential distribution. After that we estimated the power spectral density in the observed signal using the Welch 
method of spectral estimation [3^] . 

This was done in a current of InA and at the frequencies relevant for the experiment. Each spectrum was taken 
for total sampling time of 0.001 seconds. By changing the frequency we changed the average number of electrons 
for each period. For each such spectrum the power spectrum was calculated. Fig. 4 shows the result of such a 
simulation for cases relevant to the experimental situation. It is clear that the calculated linewidths are much larger 
than anticipated from lifetime considerations. The asymmetric lineshapes observed are very similar to those observed 
from the experiment (compare with Fig. 3). A rapid increase in the linewidth is predicted (Fig. 5). The data we 
have so far support this prediction but more experiments are required (Fig. 6). Using longer lifetimes gives narrower 
lineshapes but the overall behavior is the same. We want to emphasize that we do not take into account the backaction 
effect of the tunneling electrons on the precessing spin. A similar behavior is expected when the tunneling current is 
reduced. However this may affect other things, and we think that the dependence of the linewidth on the field is a 
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FIG. 4: The simulated lineshapes for a periodic function that is accumulated with a Poisson distribution of sampling times. A 
current of InA was assumed. The lineshape is clearly asymmetric. 
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FIG. 5: The predicted linewidths as a function of the magnetic field. 

more informative measurement. 

The lineshape and linewidth which are shown in Figs. 4 and 5 were calculated with the 1 /f spin noise model, but 
the same results will be observed from other models as well. Sampling of a sinusoidal function at times given by an 
exponential distribution will give similar results. 

It is clear that the calculated linewidths are at least one order of magnitude larger than the experimental linewidths. 
This difference reveals that the exponential distribution of sampling times is totally unrealistic. Because of the small 
size of the tunneling region in the STM tip, two electrons can not tunnel at the same time because of the strong 
electrostatic repulsion between them. Thus the exponential distribution has to be modified such that some temporal 
correlations will be introduced in the tunneling times. We are currently trying to get more data on the linewidths 
in different magnetic fields, currents and bias voltages, with the hope that ESR-STM lineshapes and linewidths will 
reveal details on the temporal correlations in the tunneling process. 

Conclusion 

We have proposed a 1/f noise as a mechanism of the coupling of tunneling current from the STM tip to localized 
spin S. This mechanism allows one to detect a signal from single spin even in the case when there is no dc spin 
polarization of the tunneling current. Instead we argue that the electric current will have a contribution coming from 
the coupling of spin fluctuating current to the local spin. The best way to detect a single spin in this approach is to 
perform a difference experiment where the noise spectum is taken at the local spin site and then at the non-magnetic 
site elsewhere on the surface. The difference of two noise spectra would reveal the localized spin contribution. We 
find the signal to noise ratio for this model to be on the order S/N ~ 1. We discussed the affect of random sampling 
times on the linewidths. The results explain the assymmetric lineshapes commonly observed with ESR-STM. The 
results predict a rapid increase of the linewidth with the field. The linewidths calculated are too broad which indicates 



8 




35(1 

Magnetic Held (Gauss) 



FIG. 6: The linewidths as a function of the magnetic field. The upper curve is the linewidths measured for Si Pb center |2£ 
while the lower curve for BDPA molecules 1391 



that there must be correlations between the times of the tunneling events due to Coulomb repulsion. We think that 
lincwidth measurements can provide a lot of information on the temporal nature of the tunneling process. 
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